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Abstract. The heat transport in heavy-doped n-GaAs has been investigated 
at temperatures T = 300 K and 77 K using the irradiation of the metal- 
semiconductor contact by modulated C02-laser radiation. It is shown this 
approach giving an opportunity to determine the thermo-diffusion coefficient 
X and Seebeck coefficient St without direct measurements of the temperature 
gradient. It was also found out that the thermalization length of hot electrons 
exceeds in of order of magnitude the assessment which can be done based on the 
reference data for GaAs. To elucidate the origin of the observed phenomenon the 
measurements were conducted out with Schottky contacts made on the thin doped 
GaAs layer epitaxially grown on the semi-insulating GaAs substrate. In this case 
the degenerate electron gas occupies only insignificant part of the heat-conducting 
medium. In addition, the injection of hot electrons into the semiconductor 
by current pulses through the Schottky barrier was used to clear up whether 
there is a dependence of the effect on the method of the electron heating. The 
nonequilibrium of LO-phonons and the change in the electron-phonon collisional 
integral due to the non-equilibrium pair correlations of the electrons are suggested 
as a possible explanation, f 



1. Experiment and data handling 

The frequency dependence of thermal response of the degenerate electron gas in GaAs 
(the Fermi energy ~ 80 — 130 meV) has been studied at spatially-nonuniform heating. 
The surplus energy of the electron gas was produced by free carriers absorption of 
chopped 10- fim laser radiation (see Fig|TJ. 

The opposite side of the sample was held at the equilibrium temperature. The 
change in the electron temperature T close to the metal-semiconductor interface led 
to change in the metal-semiconductor tunnel junction resistance and to the nascent 
thermoelectromotive voltage. The plot in Fig[21 demonstrates the dependence of the 
response AV on the voltage bias V measured by phase-sensitive detector. At V = 
the response is the sole thermoelectromotive signal. At V 7^ the response is the 
sum of the thermo-e.m.f. and the change in the V in the case of non-zero current I 
through the junction due to the temperature dependence of the junction resistance R. 

X Presented in part on 13th International Conference "Nonequilibrium Carrier Dynamics in 
Semiconductors" (HCIS-13), Modena, Italy, 2003 and 6th Russian Semiconductor Physics Conference, 
S-Peterburg, 2003 
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Figure 1. A draft of measurement scheme and different types of samples' layout 
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Figure 2. Bias dependence of the thermoresponse to C02-laser radiation. The 
dependence of the response voltage AV on the chopper frequency F was the 
subject of the investigation. The bias voltage Vq where the response changes in 
its sign can be used to determine the Seebeck coefficient St according to the 
Eq.0. Experimental parameters: electron density N = 7.0 ■ 10 18 cm -3 , incident 
radiation power P = 250 mW, modulation frequency F = 570 Hz, environment 
temperature T = 77 K. 
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1.1. Bias dependence of thermoresponse and determination of Seebeck coefficient 



Such an experimental setup allows us to suggest a simple way of measurement 
of the Seebeck coefficient St for the degenerate electron gas in the heavy-doped 
semiconductor when the standard methods give a large spread of values (see e.g. 
PP). A simple circuit small-signal analysis with the current-driven junction gives the 
condition for the bias Vq where AV = 0: 

AV = Vof^^j -S T AT = 0. (1) 

Here AT is the temperature difference averaged over the metal contact area. Evidently, 
it is possible to determine St from Eq. without AT measuring. All other quantities 
involved in this relation characterize the junction itself and can be easily measured. 
As an example, the value of Vq = 32 mV presented in the Fig^ being combined 
with the value of dh\R/dT ~ — 1CP 3 KT 1 gives St — —30 /zV/K in good agreement 
with the value expected for the indicated electron density and temperature. In fact, 
the derivation of Eq. is based on the potentiality of the temperature- and Coulomb 
fields and on the constancy of the Coulomb potential along all metal contacts including 
ohmic ones This allows to prove the Eq.JTJ is independent of geometry of the 
current and heat flows in the samples as long as the junction resistance is much greater 
than the bulk resistance of the substrate. Therefore, both type of the structures shown 
in Fig^have been suitable for that sort of measurements. 
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Figure 3. Comparative plot of electron density dependence of Seebeck coefficient 
in n-GaAs at T = 300 K. Selected by color data are obtained with the technique 
suggested in present work. 
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The Seebeck coefficients for severai structures with various values of electron 
density measured by this technique are presented in the Table and compared with 
known data in Fig|3] 



Table. Parameters of structures and respective Seebeck coefficients 



preparation 


structure 


N, 


R300K, 


din R/dT, 


V , 


St, 


AT-, K° 


technology 


type 


10 18 cm~ 3 


KOhm 


io- 3 /k° 


mV 


AtV/K° 


(F=600 Hz) 


MOCVD 


B 


1.27 


124.0 


-6.49 


34.4 


-223 


0,03 


MBE 


B 


3.25 


0.945 


-17.8 


6.8 


-121 


0,07 


MBE 


B 


3.95 


5.590 


-11.1 


9.93 


-110 


0,06 


zone-melting 


A 


6.80 


0.004 


-7.64 


13.4 


-102 


0,11 


MBE 


B 


2.60§ 


7.480 


-6.38 


19.2 


-123 


0,03 



1.2. Modulation frequency dependence of thermoresponse and determination of 
hot- electron thermalization length 

The measured frequency dependences of the thermoresponse for three junctions at two 
ambient temperatures are shown in Fig0]by markers. The curves 1-3 were measured 
with the same junction made on the bulk-doped GaAs substrate. The curves 4 and 
5 were measured with the junctions made on the epitaxially-doped side of the semi- 
insulating GaAs substrate. The thicknesses of the substrates varied in the 0.3 -r 0.8 
mm range. The thickness of the doped epi-layer of the sample for curve 4 was about 1 
/im. The junction relating to the curve 5 was formed by surface-deposited metal film 
and 2D electron gas in the 5-doped channel placed at the 20 nm distance from the 
semiconductor surface. 

It is noticeably that curves 4 and 5 show more slow-down behavior than the 
curve 1 as a function of the modulation frequency F in the low-frequency region up 
to F~ 10 3 Hz. In fact, this dependence is closer to 1/F 1 / 2 for the curve 4 and 5 and 
to l/F for the curve 1. More quantitative comparison of the measured curves with 
the heat transport theory has been carried out on the base of the theoretical analysis 
presented in 0]. 

For the sake of simplicity of the discussion let us consider the corresponding 
expression for the response which takes into account all essential details of the 
measurements but the finite thickness of the substrate: 

AT(O) = m^ZlK r ds eM-s^ /2 ) Ms a)J ( S d) (JJ) 



TTQK 



Certainly, the real calculations have been conducted with the full more complicated 
formula. Here -P(O) is the power of the incident radiation beam, = 2ttF, r is the 
radiation reflection coefficient, a is the radius of the metal contact, b is the dispersion 
of the Gaussian distribution of the radiation power across the beam in the focal spot, 
d is the distance along the sample surface between the centre of the metal contact 
and the focal spot, k is the semiconductor heat conductivity and K is the free carrier 
absorption coefficient. In the case of heat transport in homogeneous medium the 
quantity fcrj = y i&T/x is the root of dispersion equation. At the real f2 the Re 1/kn 
plays the role of the wavelength and the Im 1 /ka is the attenuation length of the heat 
wave. 

§ Equivalent volume density of 2D electrons with A^d = 9 X 10 11 cm -2 in <5-doped layer of order of 
3.5 nm width. 
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Figure 4. The dependence of the heating signals S on the frequency of CC>2-laser 
beam chopper for several tunnel Schottky-barrier junctions with various electron 
densities N e and diameters of the metal contact 2a (curves 1, 3-5). Curve 2 was 
measured using Q-switched CO2 laser and the posterior Fourier-transform data 
processing. Markers denote experimental data. Solid curves have been calculated 
from Eq.JJJ with parameters described in the text. 



For the inhomogeneous heating under consideration the frequency dependence of 
the right-hand side of Eq.Q at d = is determined an interrelation between \ka\ 
and reminder of the scale parameters a, b, K. The case d 7^ must be considered 
separately. 

The qualitative analysis of the expression J5| as a function of f2 with account 
for the actual limitation of the integration region by the damping exponential- or 
oscillating Bessel function (or both) in the integrand shows that at factual values 
of parameters a and b of order 200 — 250 /jm and l/K of order of 20 — 50 (im the 
response must down as ~ l/y/F in the range of 10 2 Hz< F < 10 4 Hz and transforms 
to ~ 1/F law at higher frequencies. All these estimations are made at T = 300 
K and reference value of the thermodiffusion coefficient X300 = 0.26 cm 2 /s obtained 
from thermoconductivity K300 = 0.46 Wcm~ 1 K _1 and heat capacity C300 = 1-75 
Jcm -3 K -1 of GaAs given in [3] and 0- Respective values for T = 77 K are K77 = 3.2 
Wcm^K- 1 , C 77 = 0.85 Jcm^R" 1 , X77 = 3.8 cm 2 /s . 
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The observed behavior of the curve 1-3 in Fig0] is evidently in contradiction 
with the reasoning presented above and exact computations of Eq.© support this 
evidence. Since the experimental data is represented in the form of amplitude of 
oscillating response, the thin theoretical curves in Fig ^represents the absolute value 
of the Eq.J5J. The fitted solid curves show results of calculations with optimal values 
of the thermodiffusivity coefficient x an d the reference values ^ 0] of the radiation 
absorption coefficient K. 1-2: x = 0.026 cm 2 /s, K = 450 cm -1 , 3: x = 0.38 cm 2 /s, 
K = 450 cm -1 . The theoretical curves 4 and 5 were calculated at the reference values 
of x = 0.26 cm 2 /s and K = 200 cm" . These results prove that the response of 
the samples where degenerate electron gas occupies only small part of the substrate 
volume can be described with the standard value of the thermodiffusivity x- Otherwise 
(curves 1-3), the values of x were needed to decrease in 10 times. The very sharp drop 
of the curve 4 in the high-frequency tail is the result of large displacement (250/im) 
between the centers of the focused laser beam (b — 45/im) and the metal electrode. 
Such a kind of measurements allows to determine x directly disregarding uncertainties 
of other parameters involved as it will farther be explained. 

By inspecting of Fig0] one can see that exact calculations of Eq.(j2J) with the 
indicated reference values of the parameters x an d K enable to describe the measured 
frequency dependence for the junctions on the epitaxially-grown substrates only 
(curves 4 and 5). Just in this case the heat transport passes mainly through the part 
of the semiconductor substrate that is free of electrons. In the other case, presented by 
curves 1-3, the used values of x were reduced in 10 times to reach an agreement with 
experimental data under all measurement conditions. Almost the same results can 
be obtained with the reference values of x an d reduced in 10 times K . The meaning 
of this interchangeability will be discussed later. It follows that just the electronic 
constituent is responsible for the anomalous heat propagation in the semiconductor. 

To check if the observed effect depends on the way of the electron heating, the 
pulsed current heating at the reverse- and the forward polarities was realized by 
specially designed switching electronic circuit. The circuit enabled to measure the 
thermoresponse signal of order of 200 nV in 50 — 100 [is after the voltage pulse of 1 V 
magnitude was supplied, in spite of the large capacitance ~ 2 — 5 nF of the investigated 
junctions. The time-dependent pulsed response measured by the boxcar integrator was 
converted then in the frequency domain using digital Fourier transformation. 

It turn out that frequency dependences obtained at the reverse and forward 
pulsed bias are different. At the forward bias the non-equilibrium electrons were 
generated inside the metal film and the emitted phonons only were primary mechanism 
for the heat transport into the semiconductor. In that case, the calculated 
frequency dependence of the response manifested oc 1/F 1 / 2 drop-down behavior up 
to several hundred kHz corresponding to the near-surface heat generation and the 
normally expected lattice heat transport into the semiconductor. At the reverse 
bias electrons were injected across the Schottky barrier from the metal electrode into 
the semiconductor with kinetic energy about 1 eV above Fermi level and produced 
highly non-equilibrium state of the electron gas due to intensive electron-electron 
collisions. That is the Joule heat was directly scattered inside the electron gas of 
the semiconductor. The respective frequency dependence shows the main sign of the 
abnormal heat transport, e.g., the fall off more sharp then 1/F 1 / 2 beginning from 
several hundred Hz. 

In principle, the thermoresponse can fall off with the frequency more quickly 
then 1/F even under the normal heat transport conditions as it is seen from the 
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experimental and calculated curves 4 in FigQ] That case takes place for d 7^ 0, i.e., 
when there is a distance between the centers of the metal electrode and the focal spot. 
Surely, this behavior is contained in the Eq.© and can be explained on that ground. 
However, there is a simpler way to understanding that gives us also some practically 
useful means. The solution of the heat conductivity equation far away the excitation 
region has the asymptotic form like exp(ikad)/d. Taking into account the definition of 
fcn given above it is easy to obtain the dependence of the response S on the frequency 
fl at d constant in the form 

In S cc -dy/Q/2 X - (3) 

Accordingly to this expression there is the linear dependence of In S on fl 1 / 2 . 
That dependence can be applied to the determination of x without any reference to 
the junction (heat detector) characteristics and irradiation conditions excluding the 
distance d. The application of the claimed method to the curve 4 results in the value 
of x ~ 0.23 cm 2 /s in a good agreement with the adopted value 0.26 cm 2 /s. 

2. Discussion 

The observed anomalous frequency dependence of the thermal response can be related 
either to a diminution of an effective thermal diffusivity in the system of electrons + 
LO-phonons + acoustic phonons or to a real increasing of the size of the region where 
the thermalization of the nonequilibrium electron - LO-phonon subsystem takes place 
owing to the interaction with the thermal bath of acoustic phonons. The investigated 
frequency band of the heating modulation is too low-frequency for one to expect 
diminution of the thermal diffusivity owing to the mechanism examined in |H] . The 
change in \ — k/C may be attributed either change in k or in C or in both. Since the 
heat conductivity k appears in Eq. (J2J in front of the integral besides x we can use the 
absolute magnitude of the response to distinguish the contributions of k and C in the 
observed change of thermodiffusivity x- The realization of this procedure shown that 
the main role in the change of x belong to the heat capacitance, if this change does 
exist. The joint set of heat conductivity equations for the electrons, LO-phonons, 
and acoustic phonons analyzed in 0] allows supposing that the first case might be 
considered as the result from rising the contribution of the nonequilibrium electron - 
LO-phonon subsystem to the heat capacity of the semiconductor. 

The second case may be taken into consideration because the calculated curves 
consist also with the measured data 1-3 for the bulk-doped sample in Fig0] at 
thermodiffusivity X300 = 0.26 cm 2 /s and %77 = 3.8 cm 2 /s (reference values) and 
reduced radiation absorption coefficient K — 50 cm -1 instead of K — 500 cm -1 . 
The value x — 0.26 cm 2 /s has been confirmed by direct measurements of the 
thermoresponse as a function of the intercenter distance d according to Eq.@. The 
observed phenomenon can then be interpreted as a long-range tail effect in the spatial 
response of an inhomogeneous electron gas corresponding to the known long-time 
tail in the time response of spatially-homogeneous electron gas. As it is known, the 
latter is related with the nonequilibrium corrections to the pair correlation function 
of electrons. In [7] the expression for the electron energy loss rate via electron-phonon 
interaction was derived with explicitly found contribution of the nonequilibrium pair 
correlations. However, the role of this correction in the heat transport is unclear yet. 
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